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POSITRON PRODUCTION BY PULSARS 


Peter A. Sturrock and Kile B. Baker 


ABSTRACT 

Calculations based on two specific pulsar models show that in either 
case the 511 keV garama-ray line from the galactic center can be explained 
as the result of the annihilation of positrons produced by pulsars. 
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I . INTRODUCTION 


Leventhal, MacCallum and Stang (1978) have recently reported the 

detection of 511 keV positron annihilation radiation from the direction 

, -3 

of the galactic center. Their observed flux level of (1.22 Z 0.22) x 10 

“2 “ L 43 5 •• X 

photons cm s corresponds to an annihilation rate 10 * s If the 

radiation originates at the distance of the galactic center and if posi- 
tron annihilation produces a continuum from the ortho state as well as the 
511 keV line from the para state. The authors note that current pulsar 
models generate large positron fluxes and that "the large number of radio 
pulsars within our field of view could in principle be injecting an ade- 
quate number of positrons into the interstellar medium to account for the 
feature." They also note more speculative sources of the positrons such 
as evaporating primordial black holes, a massive black hole at the galac- 
tic center, and various matter-antimatter symmetric cosmologies.* 

The aim of this article is to examine in a little more detail their 
suggestion that pulsars can yield an adequate supply of positrons to 
account for the observed line. Since any more detailed estimate must 
necessarily be mode 1-de pen den t , the following estimates will be based on 
the "polar-cap light-cylinder" (PCLC) model (Sturrock, 1971a) and a more 
recent modification, the "polar-cap force-balance" (PCFB) model (Roberts 
and Sturrock, 1972a, b, 1973), 

* Steigman (1979) has kindly pointed out to me that current evidence 
weighs against the last possibility. Matter- antimatter annihilation 
would lead to comparable number fluxes of 511 keV photons (from e” - e + 
annihilation) and gamma rays of energy 5. 100 MeV (from nucleon- antinucleon 
annihilation), whereas current evidence (Leventhal et al., 1978; Fichtel 
et al., 1975) indicates that the former exceeds the latter by more than 
order magnitude. _ 


II. POSITRON FLUX IN THE PCLC MODEL 


The total number of electrons and positrons leaving each polar cap 
is given by equation (5.6) of Sturrock (1971). Hence the same formula 
gives the total number of positrons leaving both polar caps: 


J 


e+ 


-41.4 3 13/2 -9/2 

10 B R ' P 


( 2 . 1 ) 


where R (cm) is the radius of the neutron star, B (gauss) is the mean 
field strength of the surface and P (s) is the period. 

The total number of positrons emitted during the "lifetime” of the 
pulsar is given by 


00 

V = J J e+ <t)dt . < 2 * 2 > 

O 

We see from equation (2.1) that this may be reexpressed as 

00 

N e+ = J e+ (0) / (P/P 0 >" mdt <2 * 3) 

O 

where (for the PCLC model) m = 9/2 and P is the value of P at t o 0. 

o 

There are certain assumptions involved in setting the upper limit in 
the integral of equation (2.2) as « and in assuming that, in equation 
(2.1), P is the only quantity which changes in time. Pulsars "turn off" 
due either to cessation of pair production (Sturrock, 1971b) or to magnetic 
field decay (Manchester and Taylor, 1977,* p. 164). However, the timescale 

* This book is referred to subsequently as "MT." 
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CL 

for turn-off appears to be of order 10 years, whereas we will be con- 
sidering pulsars similar to the Crab pulsar, which initially have charac- 

3 

teristic ages of order 10 years. 

In terms of the familiar "braking Index" n, 

• „2-n .. 

P * P (2.4) 


which may be expressed as 


p "- 2 




(2 .5) 


where t is the initial value of the characteristic "age" t(s) defined by 


Hence 


P = P T 


-1 



1 + (n-1) t * t 
o 


1 

n-1 


( 2 . 6 ) 


(2.7) 


By substituting equation (2.7) into equation (2.3), we arrive at the 
following estimate for the total number of positrons produced by a pulsar: 


N 


m-n+ 


I V (0) T c 


( 2 . 8 ) 


For the braking index, we adopt the value n » 2,5 found by Groth (1975) 
for the Crab pulsar. (The changes in our estimates which would be made 
by adopting n = 3 would be negligible.) Hence, for the present model, 
(2.1) and (2.8) lead to 


v -41 .9 3 13/2 -9/2 

N , =10 BR'P't 

e+ o o 


(2.9) 


By calculating the magnetic torque and by estimating the torque also 
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from the characteristic age t, we may obtain (Sturrock, 1971b) the follow- 


ing expression for the magnetic field strength 


B = 10 


15 


1 I 1 / 2 If 3 P 


.- 1/2 


( 2 . 10 ) 


2 

where I(g cm ) is the moment of Inertia. We may now eliminate B from 
equation (2.9) by means of equation (2.10) to obtain Instead 


v . io 3 - 4 j 3 / 2 »- 5/2 p - 3 / 2 ,-y 2 

o o 


(2.11) 


In terms of the neutron-star model of Baym, Pethlck and Sutherland 
(1971), we may conveniently express the principal parameters of a neutron 
star as follows (Sturrock, Baker and Turk, 1976) : 


M 


in 33 * 45 

10 


t m 44 ' 79 
I = 10 4, 


„ 10 5 - 85 


( 2 . 12 ) 


where the maximum mass is given by y = 1. Then equation (2.11) takes the 
form 


l0 56.0 11/4 p -3/2 -1/2 

^ o o 


(2.13) 
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III. PCFB MODEL 


It Is possible to repeat the calculations of Sturrock (1971a) in 
terms of the PCFB model. Then equation (2.1) is replaced by 


J = icf 100 m- 5/s b 3 r 13/2 p- n/3 
e+ 


(3.1) 


Since now m = 11/3, equation (2.9) is replaced by 


„ 10.3 -5/6 3 13/2 -11/3 

N =10 M B R / P ' t 
e+ o o 


(3.2) 


On noting that, for the present model, 


2.5 1/3 1 / 2 -3 2/3 -1/2 

B = 10 M v I R P ' t 


(3.3) 


and substituting this expression in equation (3.2), we obtain the following 
expression for the total number of positrons produced by a pulsar according 
to the PCFB model: 


N 

e+ 


10 


- 2.8 


1/6 3/2 -5/2 

M ■ I ' R ' 




(3.4) 


On using the paramaterization of equation (2.12), this formula becomes 


N 


e+ 


1q 55.3 35/12 -5/3 

u o 



(3.5) 
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IV DISCUSSION 


Estimates will be based upon the Crab pulsar, since this Is the one 

~ 1 48 

for which we have most information. At the present time, P = 10 ’ and 
t = IQ 10 ’ °. However, on noting that the Crab pulsar was "born’' 925 y 

(10 ’ s) ago, we find [by using equation (2.7) to work backwards] that 

„ -1.71 . . n-1 10 . 54 

P r 10 and, since t * p , t =10 

o o 

The total power budget of a pulsar is given by 


_ .2 . -2 -1 
= ( 2 jr) IP t 


(4.1) 


which may be expressed approximately as 


o . _46 . 4 -2 -1 

S T = 10 (i P t 


(4.2) 


38 “1 

Hence the estimated power budget of the Crab nebula 10 erg s (MT, p.60) 

-0.5 

indicates that u sj 10 .We may now estimate the total number of posi- 
trons ejected by a pulsar such as the Crab by using the above estimate of 

the mass, the above estimates of the initial values of P and t, and 

51 .9 

equations (2.13) and (3.5). These estimates are found to be 10 for 

51 4 

the PCLC model and 10 for the PCFB model. In view of the uncertainty 

in our estimate of the mass (not to mention uncertainties in the models), 

this difference is insignificant. We may therefore Ignore the differences 

51.7 

between the two models from this point on and adopt the mean value 10 

for the total number of positrons produced by a pulsar such as the Crab. 

43 5 - 1 

On referring back to the annihilation rate 10 * s in the galactic 

center inferred by Leventhal et al. (1978), we see that the positrons can 
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be supplied by pulsars if the birthrate of pulsars in the galactic center 

, “8.2 -1 -0.7 -1 

is 10 s or 10 y 

Manchester and Taylor (MT, p.138) estimate the (area) density of 

"visible" pulsars in the solar neighborhood to be 90kpc or about 
-4,0 -2 

10 pc . The distribution of pulsars in the perlod-pulse-width dia- 
gram (MT, p, 18) indicates that, if the typical pulsar beam is a pencil 
beam of circular cross-section, the half-angle if of the beam is 5° or 
10 radian. Since on the polar-cap model there are two such beams and 

each beam will sweep out a solid angle of 2jf.2ij( (for the orthogonal case), 

*0 8 

we may infer that we see a fraction 10 * of all pul&ars. Hence the 

—3 2 -2 

density of pulsars is inferred to be 10 pc . The average age of 


pulsars is estimated on kinematic considerations (MT, p.162) to be about 
6 -9 2 -2 

10 y. This leads to 10 * pc y as an estimate of the birth rate of 


pulsars. 

A current estimate of the mass density in the solar neighborhood 
- 1 -3 

(Whitley, 1977) is 10 M 0 pc . The mass distribution with respect to the 

galactic plane (Allen, 1973) has an equivalent thickness of 350 pc or 

2.5 15 —2 

10 pc, so that the mass density in the solar neighborhood is 10 ' M^pc 

Hence we infer that the mean birthrate of pulsars, referred to mass rather 

than area, is 10 M 0 y 

The entrance aperture of the Leventhal et al. instrument has an FWHM 

width of 15°. Adopting this as a sharp edge to the window, we see that 

the instrument accepts all emission originating closer than 1.3 kpc from 

the galactic center. Extrapolation of the data presented by Oort (1977) 

10 3 

indicates that the total mass within this circle is 10 ' M^. Hence the 

-0 4 - 1 

pulsar birthrate in this region should be about 10 * y 
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The above estimate of the birthrate In the region of the galactic 
center viewed by the Leventhal Instrument differs only by a factor of 2 
from what Is required to explain their 511 keV gamma-ray flux as resulting 
from positrons produced by pulsars. Although the closeness of the agree- 
ment should not be taken too seriously, since the calculation of the re- 
quired birthrate Is sensitive to the assumed pulsar magnetic field strength, 
rotation rate and mass, it seems fair to conclude that, within the limi- 
tations of current knowledge, the Intensity of the 511 keV gamma-ray line 
from the galactic center can be explained as the result of annihilation 
of positrons produced by pulsars. 

The above conclusion is necessarily model-dependent. For comparison, 

we may note that the model of the Crab pulsar developed by Cheng and 

36 * 1 

Ruderman (1977) leads to a flux of secondary positrons of about 10 s 
[According to the parameters of that model, it appears that photons radi- 
ated (by the synchrotron mechanism) from the secondary particles will 
have too low energy to lead to additional stages of pair production.] 

This flux is smaller by a factor of about 10 5 than the flux produced on 
either the PCLC or PCFB model, indicating that one could not explain the 
511 keV line from the galactic center as the result of positron production 
by pulsars on the basis of the Cheng-Ruderman model. 
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